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GEOPHONE COUPLING 

The present invention relates to seismic geophone coupling, and in particular to 
geophone coupling in seismic surveys conducted at the sea floor. 

There are a number of methods that can be used when conducting seismic 
surveys at the sea floor. Generally, a vessel at the surface activates a signal source 
immersed in water, which generates a pressure wave in the water. An array of seismic 
sensors, such as a Nessie™ 4C multiwave array, or one or more Ocean Bottom Cables / 
Seismometers (OBC/OBS) is provided on the seabed. The term 4C here indicates 4 
component, because the sensors detect the reflected P-waves and the X, Y & Z 
components of the reflected shear waves. The OBC has a number of multicomponent 
receivers or receiver groups, consisting of geophones, which measure, among other 
components, the horizontal velocity of the sea floor in two directions, ^(inline with the 
cable) and Y (crossline to the cable). The signal from the geophones is then usually 
recorded on a vessel at the surface. 

The signal generated by the source initially propagates through the water as a 
longitudinal wave, known as a P-wave. This wave will propagate through the sea, and 
then through layers under the sea bed. After the firing of the source, the OBC will 
record the arrival of the '\vater break" or direct wave, followed by reflections from 
interfaces such as the water surface, the sea floor and layers under the sea floor. 
Depending on the angle of incidence, mode conversions can occur at each interface. 
Thus the energy of the wave may propagate through the material under the sea bed 
partly in the form of a longitudinal P-wave, and partly in the form of a transverse or PS- 
wave. The PS-wave is largely visible in the horizontal X and Y components measured. 

It is known that such systems can suffer from poor sensor coupling in certain 
circumstances, and different geophone response and coupling can arise for different 
components. The X-component coupling for a PS-wave of a Nessie™ 4C multiwave 
array deployed on a hard sea bed is known to be the least reliable component. 
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It is shown in Krohn, Chr., 1984, Geophone Ground Coupling, Geophysics 49, 
pp. 722-731, that poor coupling of geophones can be explained using a model for the 
geophone ground coupling. The geophone ground coupling is modelled as a damped 
oscillator. 

U.S. Patent No. 5,235,554 (Barr & Sanders) describes a method for correction of 
differences in impulse response between the Z-component geophone and a hydrophone 
using water breaks. 

U.S. Patent No. 5,724,306 (Barr) presents a correction method for the Z- 
component using hydrophone measurements and a model for geophone response. In an 
inversion procedure differences in transfer functions between the sensor and the model 
are minimized by adjusting the resonant frequency and damping parameters of the 
model. 

U.S. Patent No 6,021,090 (Gaiser, Barr and Paffenholz) presented a method for 
correction of the /-component of OBC data using the Z-component. His method 
minimises the energy on the transverse-horizontal component of first breaks and early 
near-offset arrivals. The PS-waves are later arrivals on larger offset shots. 

According to a first aspect the present invention provides a method of analysing 
a seismic signal comprising two orthogonal horizontal components recorded by a sensor 
package containing two horizontal geophones, the method comprising generating a 
frequency dependent calibration operator to correct data corresponding to one 
component using data corresponding to the other component in order to compensate for 
different coupling between the geophone and each component of the signal. 

From here on in this specification a single sensor package or sensor group with 
one output is considered to design the calibration operator to be applied to compensate 
for inconsistent coupling at a location. Extension of the invention to a real survey 
including several sensor packages is straightforward because the operators are designed 
and applied in a receiver consistent manner. 
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Preferred features of the invention are set out in the accompanying dependent 

claims. 

According to a second aspect, the invention provides a method of performing a 
seismic survey of earth formations beneath the seabed, comprising generating a signal, 
measuring the signal at the seabed using a geophone, and analysing the signal as 
described above. 

Preferred embodiments of the invention provide a means of compensating for 
inconsistent /--coupling without the need for any modelling of the behaviour of the 
geophone as a damped oscillator, and without the need for determining any correlation 
between the behaviour of the Z-component and the horizontal components. 

Some preferred embodiments of the invention will now be described by way of 
example only and with reference to the accompanying drawings, in which: 

Figure 1 shows schematically the elements of a multicomponent seismic survey at the 
sea floor, 

Figure 2 shows the geometry of various signals arriving at an Ocean Bottom Cable 
(OBC); 

Figure 3 shows the output from a well coupled sensor and a poorly coupled sensor in 
response to a signal at 45° to the ^-direction; 

Figure 4 shows the geometry of a signal arriving at a geophone at angle 0\o the X- 
direction of an OBC; 

Figure 5 is a flow chart showing an algorithm for the correction of all the /"-components 
of a Common Receiver Gather (CRG) using the calibration operator designed on the fly. 



Figure 6 shows the X and Y components recorded at two neighbour receiver 
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Figure 7 shows the signals recorded by three geophones (two in the seabed plane and 
one vertical to it) and a hydrophone all four embedded in the cable at the same location. 

Figure 8 shows the Y component before and after calibration together with the X 
component at the same receiver location. 

Figure 9 shows the receiver consistent application of several calibration operators to a 
Common Azimuth Gather. 

A detailed description of the figures follow. 

Figure 1 shows an arrangement used to perform a multicomponent seismic 
survey acquired at the sea floor. On the sea bed 1 an Ocean Bottom Cable (OBC) 2 is 
deployed. The OBC 2 has a number of multi-component receivers or receiver groups 3 
comprising geophones that each measure the horizontal velocity of the sea floor 1 in 
two directions, X and Y. The geophone signal is recorded on a vessel 4 at the surface. 
While the motion of the sea bed is recorded, another vessel 5 fires a source 6, for 
example an airgun array, in the water. Following the firing of the source 6, the OBC 2 
will record the water break or direct wave, followed by signals generated by reflections 
from interfaces such as the water surface, the sea floor 1 and interfaces 7 between layers 
8, 9 under the sea floor 1. Depending of the angle of incidence, at each interface mode 
conversions can occur. The incidence P-wave 10 is shown in Figure 1 reflected from 
the sub sea floor interface 7 as a combination of a P-wave 1 1 and a S-wave 12 The S- 
wave 12 is detected by the geophones 2 mainly in the horizontal components. 

The source 6 used in such surveys is usually an airgun array, which is a 
compressional source, but any other source of seismic energy can be used such as a 
shear-wave source (on or under the seabed), marine vibrator or earthquake. Although 
the source 6 is shown in Figure 1 as being immersed in the water, the invention will 
work equally well for a source located at or under the sea floor. 
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Figure 2 shows a range of possible shot geometries. A signal 13 directed along 
the x-axis of the OBC 2 is known as an inline shot, and a signal 1 4 parallel to the y-axis 
is known as a crossline shot. A shot 15 at 45° to the x-axis is also shown. Following 
such a shot, identical signals for the X and Y component would be expected for a well 
coupled geophone. This is true under the assumption of an isotropic one-dimensional 
layered earth. 

If the geophone is not well coupled the signals for the X and 7-components may 
not be identical. In Figure 3 the signals from a well coupled geophone and a poorly 
coupled geophone are compared. Trace 16 is the ^-component of the signal recorded by 
a well coupled geophone. Trace 17 is the X-component of a signal recorded by a poorly 
coupled geophone. Trace 18 is the /-component of the signal recorded by the well 
coupled geophone, and trace 19 is the /-component of the signal recorded by the poorly 
coupled geophone. All of the traces show the signal varying with time. 

The 4 Vaterbreak" signal arrives first, after 0.5 seconds, and is shown at 20. This 
is the signal generated by the incoming P-wave directly from the source. Since this 
wave is propagated through the water it is well coupled on both geophones, which 
normally rest in the water on the sea bed. The P-reflection 1 1 (see Figure 1) arrives 
next, and is recorded at 21. This signal is also well coupled on both geophones, as even 
the P-reflection 1 1 arriving from the sub sea floor interface 7 will transmit most of its 
energy into the water across the interface of the sea bed 1. The PS-reflections 12 (see 
Figure 1) are shown generally at 22. Very little of the energy of the PS-reflections 12 
can be transmitted into the water so the coupling of the geophones to the sea bed is now 
crucial. 

The A-components 16, 17 of the PS-reflections 22 recorded by the two 
geophones are well in agreement. However, /-component signals 18, 19 recorded by 
the two geophones are different. The signal 19 recorded by the poorly coupled 
geophone is weaker that that 18 recorded by the well coupled geophone and has phase 
differences. Water break 20 and P-reflection 21 signals are therefore not representative 
for this kind of coupling behaviour. 
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Consider a signal S m 23 arriving at a geophone under azimuth 6 in the 
horizontal plane and recorded as the/ component ij=xy\ as shown in Figure 4. IT* 
geophone measures a signal proportional to the x and y component of the ground motmn 
G & , The frequency response /of the geophone is given by equation 1 where Cj is the 
coupling transfer function. 

Equation 1 

Cj (0 can vary for each component due to differences in design and degree of coupling. 
No explicit dependence of Q (f) on the angle of incidence has been expressed, in fact 
the incoming horizontal particle motion can always be decomposed in a component 
parallel to the cable (X) and in one orthogonal to it (Y). 

I„ absence of substantial azimuthal anisotropy and out of plane scattering effects the X 
and Y signals at 6=45° should be identical. In other words, if the two are equally well 
coupled (i.e. C x 0=C y (J)), the following identity would be valid: 

Equation 2 

It is assumed that if the geophone has non-identical coupling for , and y components, 
the signal recorded for they-component is multiplied by a transfer function T0) in order 
to obtain the same signal for both x and components. It is moreover assumed that tins 
transfer function is time invariant. 

G 45 .,(/)=n/)G 45 .. r 0 r ) 
Equation 3 

In absence of noise the transfer function would simply be: 

U) C y {f) 
Equation 4 
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Equation 3 can be written as: 

T(f) 



Equation 5 



A ro «arion matrix «* can be used to -o«*e in ft. xy-plane d. sensor package. 
Stated of the two horizontal geophones, by an angle p. 



v 7 [sin^ 



Equation 6 

The rotation matrix, which is defined by: 

-sin ^ 
cos _ 

Equation 7 

i ^ ♦ y Y component to simulate an ideal experiment with the two 
^nT^ 

Figure 4 shows the special case of a sensor package rotated of an angle of ^5°-*so 
that the azimuth of the rotated geophone and incoming signal is 0=45°: 
G (/) = cos <p ■ G ea (f) - sin <p • T(f) ■ G 0y (f) 
G (/) = sin <p ■ G 9a (f) + cos q> ■ Hf) ■ G 9%y (f) 

Equation 8 

to Equadon 8 the y-componen. geophone response has been coveted nsing the transfer 
function T(f). 

At *=45°, the rotated geophone responses should therefore be equal: 

Equation 9 

For a Common Receiver Garner (CRG) with wide azmruth covemge, several traces are 
available and the above formularion to obtain the trarcfer toncdon T(f> of the 
calibration fdter can be extended to a„ these traces, AT, The over-de— system of 
linear equations can be wnnen for each frequency as: 
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Equation 10 

where ^ cos^+s^ = n(a) 

Equation 11 

Emotion 1 1 also defines the Y and X arimtdhal cordon renns, which are respective* 
^and »«• * order «o have the system of equations 10 defined when Danishes, 
the system can be rewritten as: 

Equation 12 

whose least squares solution is: 

Equation 13 

For sake of notations dte above formulation to derive the calibration opendor 
has been earned in te Fonrier domain, however equation 13 expresses tha. me 
cafihration operator is dte matching fiUer becween me function n M and me 

function G*M> «* - ™* °< *" * " " " " 

rewritten in the original time domain: 

T(Z) = — * 

£G,,.,(Z)G,,., (1/Z)n 2 (^.) 

i=l 

Equation 14 

The numerator of equation 14 is the sum of rhe crosscorrelations of the X and Y 
components azimuthally co.ec.ed using respective* with the factors *„ ) and > 
The denominator is the sum of the arimumally corrected autocorrelations of the Y 
components. For efficiency reasons the denvation of the calibration operator is earned 
in the time domain. 
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Figure 5 shows the data flow for correcting the Y components using the 
algorithm described before. It is assumed that the orientation of the horizontal 
geophones has been assessed using the direct arrivals or positioning information, which 
are P waves, and are therefore less sensitive to inconsistent coupling as shown in Figure 
6 From a CRG with wide azimuth coverage 24 the water layer reverberations and the 
other P multiples are removed during a pre-processing phase. Early converted wave 
(PS) events are selected in the short to medium offset range (typically 600 to 1000m). 
The water break 20 and other early arriving signals 21 are not selected. Later arriving 
Scholte waves and mud roll are also removed. The windowed signal 26 now contains 
mainly PS-reflection energy 22. 

Next the azimuthal correction terms 27 and 28 are applied to the X and Y components. 
Finally the calibration operator, T(f), is derived using Equation 14. The corrected Y- 
component signal for each trace of the CRG is obtained by convolving the calibration 
operator with the original CRG Y traces. 

Figure 6 shows the X and Y components recorded at two neighbour locations, which 
were only 25 meters apart, labelled in the figure as Receiver station 827 and 829. 31 and 
32 are respectively the X and Y components at receiver location 827, 33 and 34 are 
respectively the X and Y components at the receiver location 829. The same plotting 
scale has been used for all these traces. 31 and 33 have comparable quality, but the 34 is 
of poorer quality than 32, the reflected signals are in fact very weak. Despite the general 
poorer quality of 34, the first arrivals 35, both direct and refracted, which are essentiaUy , 
concessional waves, have been properly recorded at 34 as well. Coupling for crossline 
geophones is typically more critical because of the smaller extent of the sensor package 
in that direction. The extreme case shown in this example highlights the need to use 
converted wave events to calibrate horizontal geophones, which is one of the claims of 
this invention. 

The effectiveness of the described calibration strategy depends on the validity of 
the assumption that inconsistent Y coupling can be compensated using only the X 
component. F.gure 7 qualitatively shows the validity of this assumption for the seabed 
data recorded with the currently available generation of seabed acquisition systems. 
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Figure 7 shows the data recorded by two horizontal (38 and 39) and one vertical 
geophone (37) embedded in a cable together with a hydrophone (36). All these 
components are assembled in the same sensor package. Because of the very low P and 
particularly S velocities of the shallow layers the incident angle, for offsets and target 
depths typical of exploration geophysics, is approximately perpendicular to the seabed 
plane. The moveout velocities of the reflections recorded by the two horizontal 
geophones should therefore be substantially smaller than those recorded by the vertical 
geophone if nocross-talk between horizontal and vertical geophones occur. Figure 7 
shows that for a typical receiver location this situation is verified. In the case of the 
cross-talk phenomenon is not negligible a hardware solution consists in assembling the 
horizontal geophones in a package separated from the vertical one. 

Figure 8 shows the horizontal components of a common receiver gather before 
(40) and after (41) calibration of the Y component. The data used to design the 
operators are black framed. 40 is the original Y. The Y mid trace has very little energy 
because has been obtained with a shot at the crosspoint between receiver and shot line, 
that is shooting on the inline. The calibration of this common receiver gather affects the 
amplitudes and phases of the Y gather. The Y amplitudes are generally scaled up and 
more comparable with the X ones (42). The X signal, as expected, substantially 
decreases at large offsets because of shooting on the crossline. 

Figure 9 shows the result of the application of the algorithm subject of this 
invention to an entire seabed seismic survey, only the traces whose azimuth is 
approximately 45 degrees are shown. Assuming that out of plane scattering effects and , 
azimuthal anisotropy have negligible effects the X and Y common azimuth gathers 
should be comparable. This is not the case with the original data, left panel 43 for the X 
and middle panel 44 for the Y. After calibration the Y traces 45 shown in the right panel 
have a frequency content comparable with the X and the resonant phenomena have been 
attenuated. 
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CLAIMS: 

1. A method of analysing a seismic signal comprising two orthogonal horizontal 
components recorded by a geophone, the method comprising generating a correction 
factor to correct data corresponding to one component using data corresponding to the 
other component in order to compensate for different coupling between the geophone 
and each component of the signal. 

2. A method as claimed in claim 1, wherein more than one seismic signal is 
measured, the method comprising using the same correction factor to correct the data 
corresponding to said one component of each signal. 

3. A method as claimed in claim 1 or 2, wherein the correction factor is determined 
using the fact that the data corresponding to the two components would be expected to 
be equal when the direction of each component is 45° to the direction of propagation of 
the signal. 

4. A method as claimed in any preceding claim, wherein the signal comprises a 
transverse PS-wave component, and wherein the correction factor is determined from 
data corresponding to the PS-wave component of the signal. 

5. A method as claimed in any preceding claim, wherein the direction of one 
horizontal component of the signal is defined as the ^-direction, this component being 
the x-component, and the direction of the other horizontal component of the signal is 
defined as the y-direction, this component being the y-component, the signal arriving at 
a horizontal angle of 9 to the ^-component, and wherein the data corresponding to they- 
component is corrected using the data corresponding to the ^-component. 

6. A method as claimed in any preceding claim, wherein the signal comprises a 
waterbreak and the direction of propagation of the signal is determined using 
polarisation analysis of data corresponding to the waterbreak. 
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7. A method as claimed in any preceding claim, wherein the horizontal angle 
between the direction of travel of the signal and one of the horizontal components of the 
signal is 9, and wherein a Fourier transform is performed on the data corresponding to 
each component of the signal, to generate a function G e . x (f) from the data corresponding 
to the x-component and a function G^(f) from the data corresponding to the y- 
component, and wherein a transfer function T(f) is generated wherein 
T(f) = ian0G^(f)/G ey (f) , the transfer function T(J) being the correction factor. 

8. A method as claimed in 7, wherein more than one signal arrives at the geophone, 
at one or more angles 9, and a Fourier transform is performed on the data corresponding 
to each component of each signal as described in claim 8, but wherein the transfer 
function T(j) is generated for the first signal only and used to correct the data 
corresponding to the ^-components of all of the signals. 

9. A method as claimed in any of claims 1 to 6, wherein more than one signal 
arrives at the geophone, at one or more angles 9 to the x-direction, and wherein a single 
transfer function is generated by which the Fourier transform of the data corresponding 
to the proponent for each signal can be multiplied in order to correct that data. 

10. A method as claimed in claim 9, wherein the transfer function is generated from 
the data from a single signal. 

11. A method as claimed in claim 9, wherein the transfer function is generated from 
the sum of data from all of the signals, the data having first been rotated through an 
angle of<p = 45° - 9. 

12. A method as claimed in claim 9, wherein the transfer function is generated from 
data from all of the signals using singular value decomposition. 

13. A method as claimed in any preceding claim, wherein the geophone is part of an 
Ocean Bottom Cable (OBC). 
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14. A method as claimed in claim 13, wherein the ^-direction is defined as being in 
the direction of the OBC. 

15. A method as claimed in any preceding claim, wherein said geophone is a sensor 
package containing two horizontal geophones. 

16. A method as claimed in any preceding claim, wherein said correction factor is a 
frequency dependent calibration operator. 

17 A method of performing a seismic survey of earth formations beneath the 
seabed, comprising generating a signal, measuring the signal at the seabed using a 
geophone, and analysing the signal using the method of any preceding claim. 



18. A 
array. 



method as claimed in claim 17, wherein the signal is generated by an airgun 



19. A method of measuring seismic data as herein described with reference to the 
accompanying drawings. 

20. A method of performing a seismic survey as herein described with reference to 
the accompanying drawings. 
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